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CHAPTER  I 


INTRODUCTION 

The  frequencies  of  signals  present  in  radio  communication  systems 
often  differ  from  design  values  as  a result  of  doppler  and/or  oscillator 
instabilities  by  amounts  which,  if  not  compensated  for,  would  preclude 
effective  signal  detection.  In  digital  communication  systems  wherein 
received  signals  are  coherently  detected,  phase  lock  loops  are  generally 
used  to  maintain  the  phases  of  signals  generated  locally  within  the 
receivers  in  alignment  with  the  phases  of  received  signals.  Differ- 
entially coherent  and  incoherent  detectors  do  not  require  phase  syn- 
chronization for  proper  operation,  but  frequency  uncertainties  must  be 
maintained  small  in  comparison  with  the  data  rate  if  such  detectors  are 
to  operate  efficiently.  Frequency  tracking  loops  are  generally  employed 
to  maintain  frequency  errors  at  acceptably  small  values  in  systems  which 
do  not  contain  coherent  detectors. 

A frequency  tracking  loop  would  ideally  operate  to  maintain  the 
difference  between  the  frequencies  of  a received  signal  (or  an  inter- 
mediate-frequency (IF)  signal  generated  by  down-converting  a received 
signal)  and  a local  oscillator  (LO)  signal  generated  within  the  loop  at 
a specified  value  by  suitably  controlling  the  frequency  of  the  LO  signal. 
Generally,  the  specified  frequency  difference  is  equal  to  either  a design 
value  for  an  intermediate  frequency  or  zero,  depending  on  the  manner  in 
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which  the  loop  is  implemented.  In  either  case,  suitable  circuits  are 
employed  to  estimate  the  amount  by  which  the  frequency  of  the  LO  signal 
should  be  changed  to  maintain  the  difference  frequency  at  the  specified 
value,  and  the  error  estimate  is  employed  to  control  the  frequency  of 
the  LO  signal  in  a suitable  closed-loop  manner.  Frequency  control  is 
generally  made  possible  by  utilizing  a voltage  controlled  oscillator 
to  generate  the  LO  signal.  In  loops  which  would  ideally  maintain  fre- 
quency differences  equal  to  design  values  for  intermediate  frequencies, 
frequency  discriminators  are  usually  employed  to  generate  voltages  which 
are  proportional  to  the  LO  frequency  errorsp].  Coherent  cycle  counting 
can  be  employed  to  estimate  LO  frequency  errors  in  loops  which  would 
ideally  maintain  frequency  differences  equal  to  zero[2,3].  Either  of 
these  approaches  can  be  utilized  effectively  in  practical  systems  -- 
particularly  in  systems  wherein  the  received  signals  have  either  constant 
or  continuous  envelopes.  However,  should  the  received  signal  have  a 
pulsed  envelope,  the  minimum  frequency  error  detectable  using  practical 
coherent  cycle  counting  circuits  would  be  approximately  equal  to  the 
reciprocal  of  the  signal-pulse  duration:  a value  which  would  be  un- 
acceptably large  in  many  applications.  Frequency  discrimination  can 
be  implemented  which  will  operate  with  reasonable  effectiveness  when 
the  received  signal  has  a pulsed  envelope,  but  such  circuits  are  analog 
devices  and  are  generally  less  than  ideally  suited  for  use  in  modern 
communication  systems. 

The  configuration  and  performance  of  a loop  which  tracks  the  fre- 
quency of  a pulsed  sinusoid  are  addressed  in  this  thesis.  A description 
of  the  loop  and  applicable  preliminary  analyses  are  presented  in  Chapter  II. 
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A linear  model  of  the  loop  is  developed  in  Chapter  III  and  an  approxi- 
mate expression  is  derived  for  the  standard  deviation  of  the  frequency 
tracking  error  due  to  the  presence  of  additive  white  Gaussian  noise  at 
the  loop's  input.  Results  obtained  by  simulating  the  loop  on  a digital 
computer  are  presented  in  Chapter  IV.  These  results  show  the  extent  to 
which  loop  nonlinearities  affect  the  standard  deviation  of  the  frequency 
tracking  error  and  the  distribution  of  the  frequency  tracking  error. 
Results  are  summarized  and  conclusions  are  presented  in  Chapter  V. 
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CHAPTER  II 


LOOP  CONFIGURATION  AND  PRELIMINARY  ANALYSES 

The  frequency  tracking  loop  investigated  in  this  thesis  is  pre- 
sumed to  be  representabl e by  the  block  diagram  given  in  Fig.  1 and  is 
designated  as  the  hybrid  frequency  tracking  loop  (HFTL).  This  designa- 
tion is  considered  to  be  appropriate  since  a loop  configured  as  shown 
in  Fig.  1 would  consist  of  analog  circuits,  digital  circuits,  and  cir- 
cuits which  would  be  partially  analog  and  partially  digital  in  nature. 

The  HFTL  is  structured  so  that  the  difference  between  the  carrier  fre- 
quency of  a desired  signal  present  at  the  input  to  the  loop  and  the  fre- 
quency of  a local  oscillator  (LO)  signal  is  estimated  periodically,  and 
the  frequency  difference  (error)  estimates  are  processed  to  generate  digi- 
tal control  words  which  are  applied  in  sequence  to  the  control  input  of  a 
digitally-controlled  variable-frequency  oscillator  (D-C/VFO):  the  L0- 
siqnal  source.*  Ideally,  the  digital  control  word  would  be  varied  so 
that  the  LO  frequency  would  be  maintained  identically  equal  to  the  carrier 
frequency  of  the  desired  signal  at  all  times. 


♦Either  a D-C/VFO  or  a conventional  analog  voltage-controlled  oscillator 
(VCO)  can  be  employed  in  a HFTL  to  generate  a variable-frequency  LO  sig- 
nal . The  preferred  approach  to  effecting  frequency  control  depends  on 
the  application  of  interest.  It  is  presumed  that  a D-C/VFO  is  employed 
to  implement  the  HFTL  for  explicitness  of  discussion. 
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Unlike  in  conventional  frequency-tracking  loops,  the  signals  pre- 
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sent  at  the  input  to  the  HFTL  is  processed  by  analog  circuits  during 
successive  intervals  in  a manner  whereby  the  input  signal  during  each 
processing  interval  is  "mapped"  into  two  sample  voltages  which  can  be 
interpreted  as  the  components  of  a two-dimensional  vector.  Generally, 
each  "vector"  can  be  considered  to  be  the  sum  of  a desired-signal  vector 
and  a noise  vector.  As  is  shown  subsequently,  the  angle  between  desired- 
signal  vectors  associated  with  two  contiguous  processing  intervals  is 
proportional  to  the  difference  between  the  desired  signal  and  local 
oscillator  signal  frequencies.  That  angle  is  estimated  by  suitably 
processing  four  composite  (signal  plus  noise)  sample  voltages.  The 
required  processing  operations  can  be  performed  most  readily  by  a digital 
processor  (computer)  as  indicated  in  Fig.  1 since  at  least  four  multi- 
plications, two  additions,  and  a division  must  be  performed  to  generate 
each  angle  estimate.  In  general,  the  angle  estimates  would  be  filtered 
to  generate  the  ditigal  control  words  which  are  applied  to  the  control 
input  of  the  D-C/VFO  in  time  sequence.  However,  in  this  thesis,  each 
"new"  control  word  is  presumed  to  be  generated  simply  by  multiplying  an 
angle  estimate  by  a constant  and  adding  the  result  to  the  preceeding 
("old")  control  word.  Loops  implemented  in  this  manner  are  compatible 
with  many  applications  of  interest. 

As  should  be  evident  from  the  foregoing  discussions,  a properly- 
implemented  HFTL  can  be  employed  to  track  the  frequency  of  a pulsed 
sinusoid  provided  the  loop  is  operated  in  conjunction  with  a subsystem 
which  generates  appropriate  timing  waveforms.  Clearly,  the  frequency 
of  a sinusoid  having  a continuous  envelope  could  also  be  tracked  by  a 
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HFTL.  However,  HFTLs  are  envisioned  as  being  used  primarily  in  appli- 
cations of  the  preceding  type;  it  is  subsequently  presumed  that  a pulsed 
sinusoid  is  to  be  tracked  by  the  HFTL.  A mathenatical  basis  for  pro- 
ceeding with  the  development  of  a linear  model  of  a HFTL  and  for  employ- 
ing a digital  computer  to  simulate  a HFTL  when  a pulsed  sinusoid  and 
noise  are  present  simultaneously  at  the  loop's  input  is  presented  in  the 
remainder  of  this  chapter. 

The  performance  of  a HFTL  depends  to  a considerable  extent  on  the 
properties  of  the  signal  and  noise  present  at  its  input.  In  this  thesis, 
the  composite  signal  present  in  the  input  to  the  HFTL  is  presumed  to  be 
expressible  as* 

(1)  Zr(t)  A Sr(t)  + nr(t) 

where  Sr(t)  represents  the  "desired"  signal  and  nr(t)  is  considered  to 
be  a sample  function  from  a stationary  stochastic  process  which  is 
Gaussian  distributed  with  zero  mean  and  white  with  a single-sided  power 
spectral  density  of  Nq  watts  per  Hertz.  Furthermore,  noise  process  n^(t) 
is  assumed  to  be  ergodic.**  The  desired  signal  is  assumed  to  be  express- 
ible as 


* The  symbol  - is  used  to  denote  equal  by  definition. 

**Occasional ly , a single  notation  is  used  herein  to  represent  both  a 
sample  function  from  a stochastic  process  and  the  family  of  functions 
which  constitute  that  process;  the  appropriate  interpretation  of  a given 
notation  should  be  clear  from  the  context  in  which  it  appears.  Here- 
after, it  is  presumed  without  elaboration  that  selected  ensemble  aver- 
ages to  be  calculated  are  each  equal  to  a time  average  which  describes, 
in  part,  the  performance  of  a single  HFTL,  i.e.,  that  selected  random 
responses  of  the  HFTL  are  ergodic  in  an  appropriate  sense  in  addition 
to  n (t)  being  ergodic. 
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where  u>  is  the  radian  frequency  to  be  tracked  by  the  HFTL , a Is  a phase 
angle  which  is  constant  at  any  value  in  the  interval  [-n,n],  and  $(t) 
represents  a (any)  phase  modulation  which  is  reproducible  within  the 
receiver.*  Signal  amplitude  A(Jt)  is  assumed  to  be  expressible  as 

oo 

(3)  A (t)  $ /2P~  l PT  (t-m  T ) 

m=-“>  hp 


where 

(4) 


and 


pT(ti 


<5>  Tr  > Tp  ' 2Ihp 

That  is,  Ar(t)  is  presumed  to  be  a train  of  nonoverlapping  rectangular 
pulses  which  are  generated  at  a uniform  rate  of  one  pulse  every  Tr 
seconds.  As  illustrated  in  Fig.  2,  each  pulse  in  A^U)  has  a Tp  second 
duration  and  an  amplitude  equal  to  It  can  easily  be  shown  that 

Pf,  is  equal  to  the  average  power  contained  in  the  desired  signal  during 
the  pulse  "on"  times. 


*For  example,  in  a digital  biphase  spread-spectrum  communication  system 
application,  t(t)  would  be  equal  to  C(t)n  where  C(t)  represents  a known 
binary  pseudo-noise  (PN)  sequence. 


Fig.  2.  Signal  amplitude  versus  time. 

The  LO  signals  present  within  the  HFTL  are  assumed  to  be  express- 
ible as 

(6)  s<x(t)  ^ Ak,  sin[u.ft  + *(t+.  ) ♦ «] 
and 

(7)  S (t)  ' Ak  + ^ ( t+*  ) + H] 

where  e represents  a timing  error  associated  with  a receive  clock,* 
is  a constant  phase  angle,  and  u>  represents  the  radian  frequency  of  the 
LO  signals.  It  is  assumed  that  all  circuits  have  unlimited  dynamic 
ranges.  One  consequence  of  this  assumption  is  that  loop  performance  is 
independent  of  LO-signal  amplitude  A^,  (provided,  of  course,  that  A,  is 
finite  and  nonzero).  Subsequently,  it  is  presumed  that  A equals  i'2/T^ 
for  convenience  of  presentation. 


*As  noted  previously,  the  HFTL  operates  in  conjunction  with  a timing 
subsystem  which  generates  all  required  timing  waveforms. 
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Input  signal  Zr(t)  is  multiplied  by  (mixed  with)  the  two  local 
oscillator  signals  in  separate  mixers,  the  product  signals  are  integrated 
over  appropriate  intervals,  and  the  integrator  output  voltages  are  sam- 
pled at  appropriate  instants  by  sample  and  hold  circuits  (see  Fig.  1). 

If  the  HFTL  is  to  operate  efficiently,  timing  error  c must  be  small  in 
comparison  with  the  reciprocal  of  the  spectral  width  of  $(t)  if  the 
spectrum  of  ( t ) is  sufficiently  "spread"  due  to  *(t)  being  nonzero. 
Should  4>(t)  be  zero,  t would  only  need  to  be  maintained  small  in  com- 
parison with  the  duration  of  the  integration  intervals.  Irrespective 
of  the  nature  of  t(t),  the  dump  time  of  the  integrators,  and  the 
aperture  time  of  the  sample  and  hold  circuits,  must  each  be  small  in 
comparison  with  the  duration  of  the  integration  intervals  if  the  HFTL 
is  to  operate  efficiently.  Throughout  the  remainder  of  this  thesis, 
parameters  e,  and  6?  are  each  considered  to  be  equal  to  zero. 

Four  sample  voltages  are  generated  for  each  pulse  in  sr(t)  as 
indicated  on  the  timing  diagram  shown  in  Fig.  3.  The  sample  voltages 
generated  by  processing  the  mth  pulse  are  designated  as  Zgx,  Zgy,  Z^, 
and  l*  These  voltages  are  defined  as  follows: 


(8) 


Z - Z (m  T + n T.  ) 
nx  x'  r hp 


and 


(9) 


7m  a 

"ny  y 


Z:..  A Z..(m  Tr  ♦ n Thp) 


where  Z (t)  and  Z (t)  represent  the  sample  and  hold  circuit  output  volt- 
x y 

ages  (see  Fig.  1).  Clearly, 
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Fig.  3.  Timing  diagram 


(10) 


,m  m . m 

Z_  = s + n 
nx  nx  nx 


and 


(ID 

where 

(12) 


,m  m . m 
Z * s + n 


ny  ny 


nx 


m a 
5 - 
nx 


hp 


mT_+nT__ 

, 4 sr(t)  s (t)  dt 

mv(n-mhpr  ex 


(13) 


sm 

ny 


'"VnThp 

"V<"-'>Thp 


sr<‘>  st,«>  dt 


(14) 

and 

(15) 


m a 
n - 
nx 


n 


m a 
j - 

ny 


"VnThp 
mV(n-,)Thp 


fmT  +nT 


nr(')  s,„<t)  dt 


r , nr(t)  se,(t>  dt 

mTr*(n"l )Thp  ^ 


It  can  be  shown  in  a straightforward  manner  that  rv  , n!"  , n?  , and  n? 

ux  Oy  lx  1 ; 

are  sample  sequences  from  independent  Gaussian  random  sequences  which 
each  have  a mean  value  of  zero  and  a variance  of  Nq/2[4], 

An  expression  for  sqx  is  obtained  by  substituting  appropriate 
expressions  for  sp(t)  and  s£x(t)  into  Eq.  (12)  and  setting  n equal  to 
zero: 
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mT  _T^  ^ sin[«rt+*(t)+a]/?7Tj^  s1n(Wtt+*(t)+e]dt 


fmT 

= </T7T~  I r (COs[(u>  -u>  )t+a-ti] 


r hP  JmT  -T. 

r hp 


- COS[((Or+u)e)t  + 24>(t)+a+t0>dt  . 


When  u>r  » n/T^p  or  the  HFTL  is  otherwise  implemented  so  that  the  upper 
sideband  signals  generated  by  the  mixers  negligibly  affects  loop  per- 
formance, it  can  be  shown  from  Eq.  (16)  that 

_ >£.  mT 

(17)  S0x  * Sirr2  Sin(Awt+a-6)  r 

hp  niTr"Thp 

where  represents  the  energy  contained  in  one-half  of  a signal  pulse, 
i.e. , 

<'8>  Ehp  4 Pr  Thp 

and  Aw  is  defined  as 


(19)  Aw  - i»>r-we 


Evaluating  Eq.  (17)  gives  the  result 


(20)  s"0x  - S cos  r* 


where 

. sin(AwT.  /2) 

<2’>  T^Tfr- 


(22)  r[J  A Aw[mTr-(T.  /2)]  +«-e 
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(23) 


(24) 


and 


(25) 


Similarly,  it  can  be  shown  that 


m c . m 
s0y  = S sin  r0 


m c m 
s1x  = S cos  rl 


m _ m 
s,  = S sin  r, 
ly  1 


where 


(26) 


J . J * , T 

rl  ' r0  AalThp 


It  is  clear  from  Eqs.  (20)  and  (23)  that  s™  and  s™  can  be  inter- 

Ox  Oy 

preted  as  the  components  of  a two-dimensional  vector.  That  is,  Eqs.  (20) 
and  (23)  are  consistent  with  defining  a signal  "vector"  for  each  value 
of  m as 


(27) 


cm  ii  m m 
-0  ' s0x  b0y  iy 


where  a and  a represent  orthogonal  unit  vectors  which  are  aligned  with 
* ~~y 

the  x and  y axes  of  a right-handed  rectangular  coordinate  frame,  respec- 
tively. Similarly,  S™  is  defined  as 

(28) 

These  vectors,  i.e.,  S^'  and  S'j’,  each  have  a magnitude  equal  to  L where 
(29)  L - |S|  - »Et 


Pm  a m . m 
S,  - s,  a + s,  a 
-1  lx  -x  ly  -y 


"hp 


sin(AioThp/2) 
Aut.  „/  2 


hp 
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It  Is  obvious  from  Eqs.  (20),  (21),  and  (23)  that  the  angle  between  sJJ 
and  unit  vector  a*  is  given  by* 


(30) 


* - 


; S > 0 
r«  + it  ; S < 0 


and,  from  Eqs.  (21),  (24),  and  (25),  that 


(31) 


r 

r1^1  ; S > 0 

r?  + tt  ; S < 0 
. I 


The  angle  between  vectors  SjJ1  and  is  designated  by  rm  and  is  given  by 


(32) 


,.m  A 


/[sf.sj]  = r1;1  - r;  = Au.  T 


hp 


irrespective  of  the  sign  of  S.  Except  for  an  explicable  sign  difference, 
this  result  was  obtained  by  Miller[5]  in  the  process  of  determining  the 
effect  of  a frequency  offset  on  the  bit  error  probability  performance 
of  a differential  detector. 

When  Thp  is  fixed  at  any  finite  positive  value,  Eq.  (32)  assigns 
a unique  value  to  the  angle  between  the  signal  vectors  for  any  value  of 
Am.  However,  if  vectors  and  S™  are  given,  the  angle  between  them 

can  be  assigned  an  explicit  value  only  if  a range  of  values  for  that 
angle  which  spans  2tt  radians  is  prespecified.  In  this  thesis,  angles 


♦The  notation  /.[X..Y]  is  employed  to  designate  the  angle  by  which  vector 
Y would  have  to  be  rotated  in  a counter  clockwise  direction  to  effect 
alignment  with  vector  X. 
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which  cannot  be  assigned  explicit  values  unless  a range  of  values  is 
prespecified  are  assigned  values  in  the  range  from  minus  u to  plus  it. 
These  values  are  designated  as  principal  values  and  are  generally  repre- 
sented by  lower-case  Greek  letters.  For  example,  the  principal  value 

of  r™,  { r1^}  , is  represented  by  y™,  i.e. 

A A pv  & 

<33>  ^ 4 {r”»Pv  4 r”  -k  • 2" 

where  k is  an  integer  or  zero  chosen  so  that  -ir  < y™  <_  u.  Since  r™ 
equals  Aw  T|)p,  angle  y™  equals  zero  when  Aw  equals  zero,  the  desired 
value,  or  when  Aw  equals  k • 2n/Thp  where  k is  an  integer  (not  zero). 
Thus,  if  estimates  of  angle  y™  for  successive  values  of  m are  employed 
as  a basis  for  correcting  the  local  oscillator  frequency  as  in  HFTLs  of 
the  type  considered  herein,  then  the  frequency-error  estimate  can  equal 
zero  when  frequency  offset  Aw  equals  k • 2ir/T^p - This  condition  is 
referred  to  as  a false  lock  and  can  occur  if  the  initial  frequency  off- 
set is  sufficiently  large  or  as  a result  of  noise  being  present  at  the 
input  to  the  loop.  From  Eq.  (29),  it  is  seen  that  the  amplitudes  of 
the  signal  vectors  equal  zero  when  Aw  equals  k • 2tt/T^p.  This  fact  can 
be  employed  as  a basis  for  identifying  and  correcting  false  locks. 

Unless  otherwise  noted,  it  is  subsequently  presumed  that  the  HFTL  oper- 
ates to  maintain  Aw  at  a value  much  smaller  than  2ff/T^p  most  of  the  time. 
Now,  define  composite  signal  vectors  Zjj  and  Z™  as 

(34)  £ = 5x  * C % * " = 0-’ 

The  magnitudes  of  these  vectors  are  represented  as 
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(35)  Zj  A |ZJJ|  ; n * 0,1 

It  is  easily  shown  that  the  composite  signal  vectors  can  be  represented 
as 

(36)  Zjj  = + n'("x  ax  + n™y  ay  ; subscript  n = 0,1 

The  angles  between  vectors  Zjj  and  associated  signal  vectors  S”1  are 
designated  by  i .e. , 

(37)  4 4g.Sj]  ; n » 0.1  , 

t 

and  are  defined  so  that  -n  < ^ <_  n.  The  angle  between  Zj]'  and  is 
define^  as  e01  and  is  given  by 

(38)  d"  ^ Zlzf.zJ]  « rj  + 
where 

(39)  y J ^ 

That  is, 

(40)  o'"  3 4u»  Thp+ 

The  principal  value  of  0m  is  designated  as  em: 

(41)  em  A {0m}pv  = om  - k • 2* 

where  k is  an  integer  or  zero  chosen  so  that  -tt  < em  n. 

A visual  representation  of  several  vectors  and  angles  which  have  been 
defined  to  this  point  is  given  in  Fig.  4.  In  this  figure,  angles  >q 


1 


Fig.  4.  Received  signal  vectors  in  the 
x,y  coordinate  frame. 

and  y™  represent  principal  values  of  angles  and  r™,  respectively, 
and  are  illustrated  in  a manner  which  is  appropriate  when  S is  positive 
(see  Eqs.  (30)  and  (31 ) ) • 

Sample  voltages  Zqx>  Zq^,  Z^,  and  Z^  are  digitized  and  input  to 
a digital  processor  in  time  sequence  (see  Figs.  1 and  3).  As  is  shown 
in  Chapter  IV,  each  set  of  (digitized)  sampled  voltages  can  be  processed 
in  either  of  several  different  ways  to  generate  a suitable  estimate  of 
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radian  frequency  error  Au>.*  However,  for  each  of  the  processing  algor- 
ithms considered,  the  frequency  error  estimate  is  a function  of  angle 
em.  When  the  HFTL  is  operating  in  the  intended  manner,  the  magnitude 
of  angle  om  is  less  than  n almost  all  the  time  and,  thus, 

(42)  am  = om  = Aw  Thp  + ^ 

almost  all  the  time.  Clearly,  the  accuracy,  with  which  Au  can  be  esti- 
mated is  dependent  on  the  properties  of  i j"’  - t1’1 . 

Expressions  for  conditional  probability  density  functions  associ- 
ated with  ij/g  and  Vjj1  are  now  derived.  To  this  end,  consider  Sg  to  be 

given  and  assume  for  explicitness  that  S is  positive.**  Introduce  a 

second  rectangular  coordinate  frame  having  axes  x‘  and  y'  oriented  so 
that  axis  x'  is  aligned  with  vector  as  shown  in  Fig.  5.  The  projec- 
tions of  vector  Zg  on  the  x1  and  y*  axes,  respectively,  are  given  by 

(43)  Zgx,  - L ♦ ngx. 

and 

<««>  Zoy‘  - %' 

where 

,.r\  m / m.  m . . m.  m 

(45)  n0x,  = (cos  y0)n0x  + (siny  0)nQy 


* It  is  assumed  that  the  quantization  error  introduced  by  the  digital 
subsystem  is  negligibly  small. 


X 

Fig.  5.  Received  signal  vectors  in  the  x',y' 
coordinate  frame. 


t 


and 


(46) 


V 


/ . rtu  m , , m.  m 
- (sin  Y0)nQ  + (cos  yQ)n0 


'0'  Oy 


Further,  the  projections  of  2^  on  the  x'  and  y'  axes,  respectively,  are 
given  by 


(47) 


’lx* 


= L cos(Aa)  Thp)  + n^. 


'lx' 


and 
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(48) 


ZTv-  = L cos(^  Thn)  + !»?„, 


where 


m / m,  m , , . m.  m 
nlx.  =(cos  Y0)nlx  + (sin  yq)^ 


(50)  nfy.  = -(sin  yJKx  + (cos  Yq)^ 

As  previously  noted,  n^,  , and  n^  are  sample  sequences  from 

independent  zero-mean  Gaussian  sequences  which  each  have  a variance  of 
N0/2.  For  this  case,  it  can  be  shown  that  njrjx,  njrjy,  n™x>  and  n™y  are 
also  sample  sequences  from  independent  zero-mean  Gaussian  sequences 
which  each  have  a variance  of  NQ/2[6].  Consequently,  the  conditional 
joint  density  function  p(  Z0x'  *z0y,/,L)  is  9iven  by 


p(Z0x 1 ,Z0y,/L) = 


exp  - 


(Z?Jv.-L)2+(Zn„.)2 


Using  this  result  and  the  transformation 


ZSx>  = ZS  cos  *0 


— 7^  pj*  > m 
= Z.  sin 
0 0 


an  expression  for  the  conditional  joint  density  function  pfZ^.^^/L)  can 
be  derived.  In  turn,  this  expression  can  be  integrated  (4)  over  z™  to 
yield 
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Similarly,  it  can  be  shown  that  substituting  i;1”  for  V’g  ’n  Eq.  (54)  gives 
the  correct  expression  for  p(iJm/L).  These  results  will  be  employed  as 
a partial  basis  for  developing  an  approximate  linear  model  for  the 
HFTL  in  Chapter  III. 

The  frequency  of  the  D-C/VFO  in  the  HFTL  is  presumed  to  be  express- 
ible in  the  equivalent  form 

(55)  io e = uio  + Kv  ec 

where  ec  represents  an  analog  equivalent  of  the  digital  control  word,  Kv 
represents  the  VFO  gain  in  radians  per  second  per  volt,  and  u)Q  represents 
the  radian  frequency  of  the  VFO  output  signal  when  ec  equals  zero:  the 
"rest"  frequency.  Control  voltage  ec  is  presumed  to  be  generated  by 
incrementing  its  value  following  the  processing  of  each  received  pulse 
by  an  amount  proportional  to  the  frequency  error  estimate.  For  conveni- 
ence of  subsequent  discussion,  differential  radian  frequencies  are 
defined  as 

(56)  Sio„  = <o„  - w 

r r o 

and 
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(57) 


6u,i.  = " “o  = Kv  ec 


The  difference  between  these  differential  frequencies  is  the  same  as 
the  difference  between  the  actual  frequencies,  i.e.. 


(58) 


6wr  - 60)^  = ur  - io0  - (w^  - w0)  = 


Au 


The  development  of  exact  expressions  in  closed  form  which  character- 
ize the  effects  of  noise  on  loop  performance  is  unmanagably  difficult. 
Thus,  resort  to  approximate  analyses  and  simulations  via  digital  com- 
putation is  necessary  to  obtain  useful  results.  These  approaches  are 
pursued  in  the  following  two  chapters. 
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CHAPTER  III 


APPROXIMATE  LINEAR  MODELiNG  AND  ANALYSIS 

An  approximate  expression  for  the  standard  deviation  of  the  HFTL 
frequency  tracking  error  due  to  noise  is  derived  in  this  chapter.  First, 
results  presented  in  the  preceding  chapter  are  employed  to  develop  an 
approximate  linear  model  of  the  HFTL.  Approximate  linear  analyses  are 
then  performed  to  obtain  the  result  of  interest. 

Two  characteristics  of  the  HFTL  virtually  preclude  the  development 
of  exact  expressions  which  describe  the  effects  of  noise  on  loop  per- 
fromance:  1)  for  given  values  of  the  sample  noise  voltages,  the  values 
of  phase  noises  ipQ  and  ^ depend  on  the  value  of  L which,  in  turn, 
depends  on  the  radian  frequency  error  Aw  (see  Eq.  (29))  and  2)  the 
observable  angle  em  equals  the  angle  being  estimated,  om,  only  when  the 
magnitude  of  the  latter  angle  is  less  than  w.  In  this  chapter,  the 
obstacles  presented  by  those  characteristics  are  avoided  by  assuming 
that  the  half-pulse  energy  to  single-sided  noise  spectral  density  ratio, 
Ehp/N0,  is  sufficiently  large  so  that  angles  Aw  Thp,  \j>g,  and  ^ are  each 
much  smaller  than  one  radian  most  of  the  time  and  that  om  is  an  observ- 
able angle  irrespective  of  its  value.  When  Ehp/N0  is  sufficiently 
large,  it  can  be  shown  from  Eq.  (54)  that  angle  0g  is  very  nearly 
Gaussian  distributed  with  a mean  value  of  zero  and  a variance  of 


(59) 


Since  and  ^ are  identically  distributed,  independent,  and  (nearly) 
Russian  processes  when  Ehp/NQ  is  sufficiently  large,  process  (which 
equals  - ^)  is  also  nearly  Gaussian  distributed  with  a mean  value 
of  a zero  a variance  of  N0/Ehp: 


(60) 


. % 


A linear  model  of  the  HFTL  which  is  consistent  with  all  foregoing 
results  when  Eh>>/N0  is  sufficiently  large  is  given  in  Fig.  6.  As  shown 
in  this  figure,  the  ''processor"  which  generates  control  voltage  ec  is 
modeled  by  an  ideal  impulse  sampler  and  an  integrator  having  a gain 
equal  to  Gp.  Each  "new"  value  of  ec  is  made  available  to  an  ideal 


Fig.  6.  A linear  model  of  the  hybrid  frequency  tracking  loop. 
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zero-order  hold  circuit  which  has  a gain  equal  to  by  a "second" 
sampler.  The  sampling  instants  of  this  second  sampler  are  delayed  by 
a processing  time  Ty  with  respect  to  the  sampling  instants  of  the  "first" 
sampler. 

The  manner  in  which  phase  noise  ym  affects  the  mean  and  standard 

A 

deviation  of  radian  frequency  error  Aw  is  now  determined  for  the  case 
where  differential  radian  frequency  6wr  is  equal  to  a constant.  Clearly, 
the  mean  value  of  the  frequency  tracking  error  due  to  noise  equals  zero 
since  y™  has  a mean  value  of  zero  and  the  loop  model  is  linear  (in  a 
sampled-data  sense).  The  variance  of  the  radian  frequency  tracking 
error  due  to  noise  is  calculated  by  first  assuming  that  Aw  equals  zero 
at  time  t equal  to  zero.  For  this  case,  it  can  be  shown  that  the  radian 
frequency  error  at  time  t equal  to  mTr  is  given  by  [7] 


I 


III-  I 

(61)  Aw(mTr)  = - GpIC  l ( 1 -G )m— 1 ~k  ; m = 1 ,2. - - * 

H k=0  A 

where  G represents  the  loop  gain  defined  as 

(62)  G 4 GpK,  Thp  . 

When  m is  sufficiently  large,  Aw(mTr)  is  stationary  and  its  variance  is 
given  by 

(63)  o2  - E-  lim  Aw(mTr ) 2 1 

Aw  1 

[m+«°  J 4 

= (G  K )2  E flim  "V  ^ (1-G)m_1"k  (l-G)01"1"*  • y*  • 

K |m+~  k=0  A=0 

o m— 1 m— 1 o o i Vo 

= (GpKv)2  lim  l l (1 -G )2m-2-k— a E{yk  y^} 

K k=0  £=0 

where  E{ x } designates  the  expected  value  of  x. 
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Now,  vk  and  v*  are  clearly  independent  when  Wi,  i.e., 

(64)  E{yJfJ)  = 0 ; k*t 

By  definition 

(65)  E((tJ)2>  4 „2 

Employing  these  latter  equations  to  evaluate  Eq.  (63)  gives 


(66) 


o2  = (G  K )2  o2  11m  "j  [(1-G)2]k 

Aw  P < BH-  k=0 

A 


(G  K )2 
p v 

6T2-GT 


0 < G < 2 


It  can  be  shown  through  a joint  use  of  this  result  and  Eqs.  (60)  and 
(62)  that  oAw  is  given  by 


A dimensionless  frequency  ratio  is  now  defined  as 


(68) 


A Af  _ Aw  T 

P ' T7T^  - 57  Thp 


The  standard  deviation  of  p is  given  by 


This  expression  was  employed  to  generate  graphs  of  op  versus  E^p/NQ  for 
several  values  of  loop  gain  G.  These  graphs  and  numerical  results 
obtained  by  simulating  the  HFTL  on  a digital  computer  are  presented 
jointly  in  the  following  chapter. 
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CHAPTER  IV 


COMPUTER-IMPLEMENTED  HFTL  SIMULATIONS 
AND  NUMERICAL  RESULTS 

A digital  computer  was  employed  to  simulate  HFTLs  having  both 
unavoidable  and  deliberately-incorporated  nonlinearities,  and  numerical 
results  were  obtained  which  show  the  effects  of  the  nonlinearities 
investigated  on  the  standard  deviation  of  frequency  ratio  p (see 
Eq.  (68)).  These  results  are  presented  in  this  chapter  and  are  compared 
with  results  predicted  by  Eq.  (69).  Data  descriptive  of  distribution 
functions  associated  with  p for  selected  loop  configurations  and  para- 
meter values  are  also  given. 

The  model  of  the  loop  simulated  by  the  computer  program  developed 
is  shown  in  Fig.  7.  Note  that,  in  this  model,  differential  radian  fre- 
quency 6^r(m)  is  accumulated  rather  than  the  control  voltage  and  sample 
voltages  defined  relative  to  the  x',y'  coordinate  frame  are  generated 
and  processed  rather  than  sample  voltages  defined  relative  to  the  x,y 
coordinate  frame.  However,  the  two  models  suggested  by  this  statement 
are  mathematically  equivalent  for  the  estimator  algorithms  to  be  con- 
sidered. A listing  of  the  simulation  program  is  given  in  Appendix  A. 

At  the  beginning  of  each  program  run,  the  differential  L.O.  fre- 
quency was  set  equal  to  zero  and  the  differential  input  frequency  was 
fixed  at  an  appropriate  value.  One-hundred  iterations  of  the  program 
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Model  of  the  hybrid  frequency  tracking  loop  simulated 
on  a digital  computer. 
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loop  were  completed  to  allow  the  simulated  loop  to  attain  steady-state 
operation.  Subsequently,  the  program  loop  was  iterated  6200  times. 

This  number  was  adequate  to  allow  accurate  estimates  of  the  mean  values 
and  standard  deviations  of  parameters  of  interest  to  be  obtained  at  an 
acceptable  confidence  level  (see  Appendix  B).  The  output  data  for  each 
run  included  the  sample  mean,  sample  variance,  and  standard  deviation 
of  the  frequency  error  normalized  to  T^p.  These  were  chosen  so  as  to 
be  compatible  with  the  theoretical  development  of  Chapter  III,  specific- 
ally, the  result  expressed  in  Eq.  (69). 

As  indicated  in  Fig.  7,  gJi(e(m))  is  related  to  Z0x”  Z0y 1 ’ *x" 

and  Z*  by  a function  f JZ^,  ,Z™yi ) . Subscript  i is  assigned  different 

values  to  designate  different  functions.  As  suggested  by  the  notation, 

each  estimation  algorithm  considered  has  a form  which  results  in 

g<(o(m))  being  dependent  only  on  the  principal  value  e(m)  of  the  angle 

between  vectors  Z™  and  Z^j.  For  convenience  of  presentation,  define 

X(Zm  ,,Zm  ,)  and  Y(Zm  ,,Zm  ,)  as 
nx’  ny'  nx'  ny 


(70) 


X(Zm  ,,Zm  Z"  .zT  , - Z?  ,Z™  , = |Zj!||Z,!,| 

nx  ny  Ox'  ly ' lx  Oy  1 —0 1 1 — 1 1 


sin  e(m) 


and 


(71) 


v / ?m  7m  \ i 7IB  7m 

YUnx,,Zny,)  " L^'L 


. 7 m 71 

Ox'Mx'  %,L}y 


m 7m  , = |zj||zj|  cos  e(m) 


Note  that  X(Z"x.,Z™  .)  is  the  z component  of  the  cross  product  of  z£J 


'nx  ’ ny 

and  Z?,  and  Y(Z™  .,Zm  .)  is  the  scalar  (dot)  product  of  Z?|  and  Z?. 
—1  nx  ny  u 

Using  this  notation,  the  five  estimator  algorithms  investigated  are 
expressible  as  follows: 


L 
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(72)  9l(e)  - MC-V)  ‘Tan 


X<znx-C> 

Y<Znx'’Zny} 


where  Tan'^  represents  a four-quadrant  inverse  tangent  function. 


(73)  g2(e)  = f 2 ( z„x • » zp y ■ ) = arctan 


X f 7m  Zm  1 

^nx 1 * nv'  ‘ 

Y(Zm  Zm  ) 

'znx'  ’ ny* 1 


where  arctan  represents  a two-quadrant  inverse  tangent  function, 

X(Zm  Zm  ) 

(74)  g3(e)  - f3(zj  ,,zj’  ,)  = — = tan  e 

3 3 nx  ny  Y(z”x,  ,2^, ) 


xfzm  zm  ) 

f,r\  ^ /-, m -,m  \ _ v nx”  ny1' 

(75)  g4(e)  = f4(znx'*zny')  = 


|v(C-.z:y«): 


and 


(76) 


95(0)  = 


fc(Zm  ,,Zm  ,) 
5'  nx  ny  ‘ 


X(Zm  ,,Zm  ,) 

nx  ny' 

rx2fZm  7™  l + Y2fZm  F 
' nx1’  ny1'  lznx' ,zny' 


= sin  9 


The  functions  g-j(e)  through  g^e)  are  illustrated  in  Figs.  8 through 
12,  respectively. 

The  standard  deviation  of  frequency  ratio  p is  shown  plotted  as  a 
function  of  the  half-pulse  energy  to  single-sided  noise  spectral  density 
ratio  for  each  g^(e)  and  for  a loop  gain  equal  to  one  in  Fig.  13. 
Theoretical  results  generated  using  Eq.  (69)  are  shown  in  this  figure 
along  with  the  numerical  data  obtained.  As  expected,  the  numerical  data 
are  in  close  agreement  with  the  analytical  result  when  the  value  of 
Ehp/N0  is  sufficiently  large  for  each  g^e).  For  the  g3 ( e ) case,  the 
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I 2 3 

e 

a limits  at  ± tan  (3r/e> 

b LIMITS  AT  ± tan  { T/4  ) 
C LIMITS  AT  ± tan  (V8  ) 
d LIMITS  AT  ± tan  (*/|6> 


Fig.  11.  The  modified  tangent  estimator  with  increased  lock 
range;  four  possible  limits  are  shown. 
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2 6 <0  14  18  22  26  30 

Efcp/N0  dB 


Fig.  13.  Standard  deviation  of  the  normalized  frequency- 
tracking error  versus  half  pulse  energy  to 
single-sided  noise  spectral  density  ratio 
for  various  estimators  and  a loop  gain 
equal  to  one:  G=1 . 


numerical  data  remain  in  close  agreement  with  the  analytical  result  as 
Ehp/N0  is  reduced  from  a large  value  down  to  10  dB.  Data  for  smaller 
values  of  Ehp/N0  are  not  shown  for  this  case  since  the  loop  did  not 
remain  locked  while  the  program  loop  was  interated  6200  times  when 
Et)p/N0  was  set  equal  to  8 dB.  The  range  of  E^p/N0  values  for  which  the 
loop  remained  locked  through  6200  program  iterations  is  indicated  in 
the  figure  for  each  ge(e).  Estimators  g4(o)  with  limits  suitably 
imposed  and  95(0)  provide  the  greatest  lock  range;  thus,  the  performance 
provided  by  these  estimators  was  determined  in  detail. 

The  effects  of  setting  the  limits  on  g4(o)  at  several  different 
values  on  op  are  shown  in  Fig.  14.  In  isolation,  these  results  suggest 
that  limits  having  a small  magnitude  would  preferably  be  employed. 
However,  the  reduction  in  op  which  can  be  effected  by  setting  the  magni- 
tude of  the  limits  at  a smaller  value  is  accompanied  by  a reduction  in 
the  effective  loop  gain.  The  limits  should  be  imposed  at  approximately 
± tan(n/4)  if  the  effective  loop  gain  is  to  remain  approximately  constant 
as  E^p/N0  is  varied  over  a wide  range. 

Data  which  show  the  effects  of  varying  the  loop  gain  on  op  for 
the  g4(o)  case  with  limits  set  at  ± tan(n/4)  and  for  the  95(0)  case  are 
shown  in  Figs.  15  and  16,  respectively.  The  numerical  results  are  in 
reasonable  agreement  with  results  obtained  using  Eq.  (69) 

A subroutine  to  the  simulation  program,  when  called,  orders  the 
values  of  p and  then  calculates  the  percentage  of  the  samples  which  have 
values  smaller  than  a specified  value.  Thus,  the  distribution  of  p was 
readily  determinable.  Distribution  data  were  obtained  for  six  cases 
of  interest;  these  data  are  given  in  Figs.  17-19.  The  straight  line 


38 


Fig.  14.  Standard  deviation  of  the  normalized  frequency  error 
versus  half-pulse  energy  to  single-sided  noise 
spectral  density  ratio  for  the  modified  tangent 
estimator  and  a loop  gain  equal  to  one:  G=1 . 
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PROBABILITY  (%) 

Normalized  frequency  error  versus  cumu.ative  probability  for  the 
modified  tangent  estimator  with  limits  at  +tan(ir/4)  and  the 
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PROBABILITY  (%) 

Fig.  18.  Normalized  frequency  error  versus  cumulative  probability  for  the  modified 
tangent  estimator  with  limits  at  itan(-/4)  and  the  sine  estimator; 
(Ehn/N.)  = 6 dB  and  G = 0.25. 


0.20 


Fig.  19.  Normalized  frequency  error  versus  cumulative  probability  for  the 
modified  tangent  estimator  with  limits  at  ±tan(ir/4)  and  the 
estimator;  (Ehn/N0)  = 8 dB  and  G=1.25. 


CHAPTER  V 


SUMMARY  AND  CONCLUSIONS 


The  configuration  and  performance  of  a hybrid  frequency  tracking 
loop  which  tracks  the  frequency  of  a pulsed  sinusoid  have  been  address- 
ed in  this  thesis.  A linear  model  of  the  loop  was  developed  and  an 
approximate  expression  was  derived  for  the  standard  deviation  of  the 
frequency  tracking  error  due  to  the  presence  of  additive  white  Gaussian 
noise  at  the  loop's  input.  Results  obtained  by  simulating  the  loop  on 
a digital  computer  were  presented  which  show  the  extent  to  which  loop 
nonlinearities  affect  the  standard  deviation  of  the  frequency  tracking 
error,  the  distribution  of  the  frequency  tracking  error,  and  the  loop 
lock  range.  These  results  show  that  the  loop  operates  with  a high 
degree  of  effectiveness  and  provide  a basis  for  designing  practical 
loops  of  the  type  investigated. 
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THE  FOLLOWING  IS  A FORTRAN  LISTING  OF  THE  PROGRAM 
USED  TO  SlHULATL  ThE  HYBRID  FREQUENCY  TRACKING 
LOOP,  THE  OUT  PUT  IS  WRITTEN  INTO  TWO  FILES. 
CuPOUT  CONTAINS  T Hf  SAMPLE  MEAN  OF  THE  LOOP 
rRROR  SIGNAL  t THE  STANDARD  DEVIATION  BASED  ON  THE 
SAMPLE  VARIANCE*  AND  POINTS  ON  THE 
aSSuCIaTEO  distribution  function 
for  a given  set  of  loop  parameters.  outclp 
PROVIDES  an  INDICATION  Of  FaLSl  LOCK. 


RAUSS  IS  A RANDOM  NUMBER  GENERATOR  THAT  OUTPUTS 
INDEPENDANT  SAMPLES  FROM  A C-AUSSIAN  DISTRIBUTION 
wlTH  THE  SPECIFIES  MEAN  AND  STANDARD  DEVIATION 

IX  IS  AN  ODD  INTEKGER  USED  AS  A STARTING  VALUE 
S IS  THE  STANDARD  DEVIATION 

AMU  IS  ThE  MEAN 
V IS  THE  OUTPUT  VALUE 
SUBROUTINE  GAUSSI IX.S.AMlUV) 

A — 0 • 0 

DO  50  IZ=1,48 
IYsIX*16645 
I F ( I Y ) 5 , 6 , 6 
I Y=I Y+6360607+1 
Y=IY 

Y=Y*,1192093E-6 

IX=IY 

A=A«-Y 

V=( A-24,0 )*.5*S+AMu 

RETURN 

END 

maIN  PROGRAM 
REAL  L 

DIMENSION  OELT A ( 6200 ) 

CALL  ASSIGN  (6HCLP0UT.6H3364S  *5) 

CALL  ASSIGN  (6H0UTCLH.6H3364S  *6) 

CALL  DEASSN 
WRITE (6,101) 

FORMAT (12X, "SNR  D8 " *14X«" MEAN M*17Xt" DEVIATION") 
DO  6 J=1 ,1 

ESTABLISH  LOOP  PARAMETERS 

X IS  THE  ENERGY  TO  NOISE'  DENSITY  RATIO  IN  DB . 

tp  is  the  integration  time 

G IS  THE  LOOP  GAIN 
X = 6. 

IF<X.EQ.-1.>  GO  TO  7 

nun=o 

N0V  = 0 
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51 

52 

53 

54 

55 

56  C 

57  C 

56 

59  C 

bO 

61 

62 

63  C 

64  10 

65  11 

66 

67  6 

6b 

69  9 

70  C 

71  2 

72 

73  4 

74 

75 

76 

77 
7b 


81 

82 

83  C 

84  C 

85  C 

86  C 

87  C 
86  C 

89  C 

90  C 

91  C 

92  C 

93  C 

94  C 

95  C 

96  C 

97  C 

98 

99 
100 


EP3l0.**(X/10.  ) 

TP*5. 137UE-5 
Gsl./TP*0.25 
W (.01x0*0 

Pl=-ATAN2<0«  * "1 • ) 

THE-  FOLLOWING  STATEMENT  DETERMINES  THE  NUMBER 
OF  ITERA»IO*'S  FOR  THE  LOOP 
DO  3 K=lt6200 

WC  IS  THE  INPUT  FREQUENCY  FUNCTION 

WC=1 • OE+4 

DELTA! K)=WC-WL01 

OFST»(WC-WL01)*TP 

detect  FALSE  LOCK 

IF (OFST.LT. -PI ) 00  TO  8 

IF(OFST.GT.PI)  GO  TO  9 

GO  TO  2 

NUN=NUN+1 

GO  TO  2 

NOVsNOV+l 

FORM  SAMPLES  used  for  input  TO  the  estimate  routine 
L=SQRT(EP)*SlN(0FST/2. )/(0FST/2. ) 

S=SQRT ( .5) 

CALL  GAUSS(4629«S*0.OtRl ) 

21=L+K1 

CALL  GAUSS(4629«S»0t0,R2) 

Z2*R2 

CALL  GAUSS(4629«S»0.0tR3) 

Z5=L*C0S  (OFST>+R3 

CALL  GAUSS(4629.St0.0iR4) 

Z4*-L»SIN(0FSTH-R4 

AGlsZl*Z4-23*z2 

AG2=Z1*Z3+Z2*Z4 

THE  FOLLOWING  STATEMENT  IS  THE  ESTIMATE  FUNCTION 
THE  EQUATIONS  FOR  VARIOUS  ESTIMATORS  ARE  A s FOLLOWS 
THE  LINEAR  ESTIMATOR 
W01aATAN2(AGl,A62) 

THE  LINEAR  ESTIMATOR  WITH  DECREASED  LOCK  RANGE 
WOl  = ATA|\|(  AG1/AG2) 

THE  TANGENT  FUNCTION  ESTIMATOR 
W01sAG1/AG2 

THE  TANGLNT  FUNCTION  ESTIMATOR  WITH  INCREASED 
LOCK  RANGE 
W01sAG1/ABS< AG2> 

FOR  THE  ABOVE  ESTIMATOR  WITH  LIMITS  ON 
TnE  OUTPUT  USE  SUBROUTINE  CLAMP 
THE  SINE  FUNCTION  ESTIMATOR 
WOl*Aul/SQRT< AG1**2+Ab2**2) 

CALL  CLAMP(AG1iAG2»W01«PI  ) 

W02*-W0l 

WL02*WL01 
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101 

tol.01s*L02+G*W02 

102 

SUH2«0 . 0 

1U5 

3 

sum=o , 0 

104 

C 

CALCULATE  STATISTICS 

105 

Du  50  K=100.fe200 

10b 

50 

SUM2=SUM2*UELTA(K  ) 

107 

fANN=SUM2/bl00. 

106 

no  5 K* 1 0 0 » b 1 0 0 

109 

5 

SuM  = SUM-M0LLTA<K>-FANN)**2 

110 

t AN»SUM2/6l0  0.*I»;’/'2/Hl 

111 

VARnSUM/blOO. 

112 

niV*SGRT(  VAK  >»TP/2/Pl 

H3 

100 

FORMAT (3F20.6) 

1 Is 

WKITtlb.-)  X. NUN. 90V 

115 

b 

-KIT£(5,100)  XtEAN.UlV 

lib 

TALL  STATSl  ( OFLT  A . TP  ) 

117 

7 

CALL  £XIT 

lift 

FND 

119 

C 

CLAMP  Is  UstD  As  THE  Es T I M A TF  kOUTINE 

WHEN 

1 2 1 

c 

LIMITS  ARE  PLACtO  ON  THE  OUTPUT  OF 

the 

TANGENT 

121 

c 

FUNCTION  estimator  WITH  INCRFASFD  i 

LOCK 

RANGE 

1 22 

c 

Abl  IS  THE  MAGNITUOE  OF  THE  VECTOR 

CROSS  PROOUcT 

123 

c 

Ab2  IS  the  vector  dot  product 

124 

c 

WOl  IS  THE  output  ESTIMATE 

125 

SUBROUTINE  CLAMP( A01. AG2.w0l.PI ) 

12b 

ALIMI T=Pl/4. 

127 

TSTLIM=AG1/ABS<AG2> 

125 

TF (TSTLIM.GT.SINT ALIMI T) /COS (ALIMI T) ) 

GO  TO  50 

129 

TF  < TSTLIM.LT. -SIN ( ALIMI T ) /COS (ALIMI T> ) 

GO  TO  51 

130 

WOlsTSTLlM 

131 

GO  TO  52 

132 

50 

W01=SIN(ALIMIT)/'C0S(ALIMIT) 

133 

GO  TO  52 

134 

51 

W01=-SIN( ALIMI T) /COS ( ALIMI  T) 

135 

52 

return 

13b 

END 

137 

C 

STATSl  IS  USED  TO  CALCULATE  POINTS 

ON 

the 

138 

C 

DISTRIBUTION  FUNCTION  OF  THE  ERROR 

OATA 

139 

c 

DELTA  IS  THE  ARRAY  OF  ERROR  VAlUFS 

190 

c 

TP  IS  THE  INTEGRATION  TIME 

141 

SUBROUTINE  STaTSI  (DELTA, TP) 

14*? 

DIMENSION  DELTA(b20O) , BOX ( 20 , 2 ) 

143 

c 

SORT  DATA  BY  SIZE 

144 

4 

FLG  = 0. 

145 

DO  2 INOEXl=100.bl99 

14b 

TNPEX2=INDBX1+1 

147 

IF(UELTA< IN0EX2) ,GE.OELTA( INnEXl ) ) 

GO 

TO  2 

1 4b 

TEMPaDCLTAT INOEX1 ) 

149 

DELTA! lNOEXl>*D£LTAi JNDEX2) 

ISO 

DLLTA( IN0EX2)=TEMP 

50 


191 

FLS«FL&+1, 

192 

2 

CONTINUE 

193 

IF(FLG.NE,0. ) GO  TO  4 

154 

C 

QUANTIZE  DATA 

195 

DO  9 1=1,20 

156 

9 

POX ( I « 2 ) =0  • 

197 

INBXsl 

156 

SMs-DLLTA(lOO) 

159 

ChsAMAXKSMt  DELTA  (6200)  > 

160 

WRl TE ( 5 « 101 ) OFLI  A ( 100  ) , OELTA  < o200  ) ,SM,CH 

161 

101 

FORMAT ( 4E 13 • 4 ) 

162 

DO  1 INDLX3=100«6200 

163 

12 

EDPTs-CH-HNBX*.l*CH 

164 

BOX( INBX.1)sEDPT»TP/2,/(.aTAN2(0. *-l. ) ) 

165 

IF <DELTA< INUExS) .GT.EDPT)  GO  to  3 

166 

5 

BOX( INBX.2)»«OX(INBX,2)*l. 

167 

GO  TO  1 

16b 

3 

INBX*IN8X+1 

169 

IF ( INdX.GT.20 ) GO  TO  6 

170 

GO  TO  12 

171 

1 

CONTINUE 

172 

C 

CALCULATE  DISTRIBUTION  POINTS 

173 

6 

RUNTOTrO. 

174 

00  7 1=1,20 

175 

BQX(I,2)=B0X(I«2> ♦RUNTOT 

1 fb 

7 

RUNTOT=BUX( I ,2) 

177 

DO  8 U=1.20 

17b 

8 

B0X(U,2)=BOX(U*2) /runtot*ioo. 

179 

DO  10  U2*l,20 

160 

10 

WRITE (5, 1000)  B0X(j2,l J ,BUX( J2,2) 

161 

1000 

F0RMAT(E12.4,F10,3) 

162 

END* 

51 
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APPENDIX  B 

CALCULATION  OF  THE  NUMBER  OF  SAMPLE  VALUES 
NEEDED  FOR  ACCURATE  SIMULATION  RESULTS 


The  number  of  sample  values  of  radian  frequency  error  Au)(m)  needed 
to  obtain  an  accurate  estimate  of  the  variance  of  Aw(m)  is  determined 
in  this  appendix  through  the  use  of  confidence  interval  theory.  It  is 
assumed  that  Aw(m)  is  a Gaussian  random  variable  with  a mean  u.  and 

Agj 

p 

variance  . It  is  known  that  the  variable  defined  by 


(77) 


* 2 
l [Aui(m)  - Aw] 

v A m=  1 
*k  2 

°Au> 


where 


(78) 


k 

l Auj(m) 

— A m=l 
Auj  r 


and  k denotes  the  sample  size  is  a chi  square  variable  with  k-1  degrees 
of  freedom.  If  k is  large,  then  the  variable  defined  by 


X,  - k+1 

(79)  Yk  4 
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is  Gaussian  with  zero  mean  and  unit  variance.  An  interval  with  end 


points  a and  b is  defined  such  that  Yk  will  be  found  in  this  interval 
99%  of  the  time: 

(80)  Pr[a  < Yk  < b]  = 0.99 

Using  standard  tables,  the  values  of  a and  b are  determined  to  be 

(81)  -a  = b = 2.576 


Thus , 

(82)  Pr[-2. 576  < Yk  < 2.576]  = 0.99 

By  using  Eqs.  (77)  and  (79)  in  the  above  expression  an  interval  con- 
taining the  sample  variance  with  probability  0.99  can  be  expressed  as 


(83) 


k-1-2.576/2(k-l) 

k 


°a  < V2  < 
Au) 


p 

where  \r  is  the  sample  variance  given  by 


k-l+2.576/2(k-l ) 
k 


2 

a 

Alo 


= 0.99 


(84) 


v2  = juLL 


o 

l [Aw(m)  - aITt 


The  desired  maximum  error  in  this  sample  variance  is  now  specified  to 
be  ±5%  of  the  true  variance.  The  desired  interval  becomes 


(85)  0.95  a2  < V2  < 1.05  a? 

' ' Au) 
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Finally,  the  desired  interval  end  points  are  equated  to  the  confidence 
interval  end  points  and  the  resulting  equations  are  solved  for  k. 
Choosing  the  largest  result  gives 


(86)  k = 5348 

This  result  is  interpreted  to  mean  that  at  least  5348  samples  must  be 
used  so  that  the  sample  variance  will  not  be  in  error  by  more  than  ±5% 
for  99%  of  the  variances  calculated. 
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